APS/123-QED 



Quantum criticality in the cubic heavy-fermion system CeIn 3 _ z Sn 2 
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We report a comprehensive study of Celns-xSn^ (0.55 < x < 0.8) single crystals close to the 
antiferromagnetic (AF) quantum critical point (QCP) at x c ~ 0.67 by means of the low-temperature 
thermal expansion and Griineisen parameter. This system represents the first example for a cubic 
heavy fermion (HF) in which Tn can be suppressed continuously down to T = 0. A characteristic 
sign change of the Griineisen parameter between the AF and paramagnetic state indicates the 
accumulation of entropy close to the QCP. The observed quantum critical behavior is compatible 
with the predictions of the itinerant theory for three-dimensional critical spinfluctuations. This has 
important implications for the role of the dimensionality in HF QCPs. 

PACS numbers: 71.10.HF,71.27.+a 



Non-Fermi-liquid (NFL) properties are observed in 
many heavy-fermion (HF) systems and frequently at- 
tributed to a nearby quantum critical point (QCP) Q. 
A QCP can arise by continuously suppressing the transi- 
tion temperature Tn of an antiferromagnetic (AF) phase 
to zero, e.g. by chemical or applied pressure or an ex- 
ternal magnetic field. QCPs are of great current interest 
due to their singular ability to influence the finite tem- 
perature properties of materials. Heavy-fermion metals 
have played the key role in the study of AF QCPs. The 
essential question is how the heavy quasiparticles evolve 
if these materials are tuned from the paramagnetic into 
the AF ordered state. The traditional picture describes a 
spin-density-wave (SDW) transition and related, a mean- 
field type of quantum critical behavior. Here, the quasi- 
particles retain their itinerant character 00- Unconven- 
tional quantum criticality which qualitatively differs from 
the standard theory of the T = SDW transition, may 
arise due to a destruction of Kondo screening. Here, the 
quasiparticles break up into their components: conduc- 
tion electrons and local 4/ moments forming magnetic 
order 0, . This locally-critical picture leads to a num- 
ber of distinct properties, including stronger than log- 
arithmic mass divergence, lu/T scaling in the dynamical 
susceptibility and a large reconstruction of the Fermi sur- 
face. Such behavior has been found at least in some HF 
systems 0, 0, Q . The central question is to identify the 
crucial parameter leading to the different types of QCPs. 
Of particular importance should be the dimensionality 
of the magnetic fluctuations, which could be reduced by 
the presence of frustration, ft is proposed in Q| that 
for magnetically three-dimensional (3D) systems without 
frustration the itinerant SDW picture should apply. On 
the other hand, 2D magnetic systems should be described 
by a locally quantum critical picture |4j. However, sys- 
tems currently under investigation, are either tetragonal, 
e^ CeNi 2 Ge 2 , YbRh 2 (Sii_ x Ge a; }2 G3 and CeCu 2 Si 2 
jllj . hexagonal, e.g. YbAgGe 12] or monoclinic, e.g. 



CeCug-xAua; [f| and the lower crystallographic symme- 
try could result in fluctuations with reduced dimension- 
ality. Therefore the dimensionality of the critical spin 
fluctuations clearly needs to be substantiated by inelas- 
tic neutron scattering experiments. In order to avoid 
this constraint, experiments on cubic systems close to 
QCPs are particular interesting. Celii3- x Sn Xl with a cu- 
bic point symmetry of Ce-atoms in the CU3AU structure 
(compare the inset of Figure 1), is thus an excellent can- 
didate for such a study, as here low-dimensional spin fluc- 
tuations can be ruled out. Thus, the interesting question 
arises, whether the mechanism of NFL-behavior in this 
system can be described by an itinerant 3D SDW theory. 

In this Letter, we present thermal expansion measure- 
ments and a Griineisen ratio analysis performed on single 
crystalline samples of the cubic system Celna-^Suj; close 
to the critical concentration, x c = 0.67, where Tn is sup- 
pressed to zero by doping. Recently, it has been shown 
that the thermal volume expansion (3 = V^ 1 (dV/dT) (V: 
sample volume) is particular suited to probe quantum- 
critical behavior, since, compared to the specific heat.it 
is much more singular in the approach to the QCP |9(. 
As a consequence, the Griineisen ratio r ~ (3/C of ther- 
mal expansion, /3(T), to specific heat, C(T), is divergent 
as T goes to zero at any pressure-sensitive QCP and the 
associated critical exponent can be used to distinguish 
between the different types of QCPs. In the itinerant 
scenario the divergence T cx I /T £ is given by e = 1/vz 
9] with v, the critical exponent for the correlation length, 
£ oc \r\ v (r: distance from the QCP) and z, the dynami- 
cal critical exponent in the divergence of the correlation 
time, r c cx £ z . For a 3D AF QCP v = 1/2 and z = 2 
yielding e = 1. Thus, a study of the Griineisen ratio 
can prove the validity of the 3D SDW picture in the title 
system provided that e = 1. 

The magnetic (x, T) phase diagram of polycrystalline 
Celna-^Sna; has been widely studied for < x < 1 by 
susceptibility 0, specific-heat 0] and resistivity mea- 
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Figure 1: Magnetic phase diagram for cubic Celns-^Sria; (x < 
1). Closed circles and diamonds indicate Tn, determined from 
specific heat [lH | and electrical resisitivity [Tr| measurements, 
respectively. Open diamonds mark T* , the upper limit of 
Landau Fermi-liquid behavior, e.g. Ap(T) oc T 2 [T3 |. Open 
triangles indicate first-order transition Ti |14|. 

surements 0, (see Figure 1). Whereas Tjy for un- 
doped Celn3 vanishes discontinuously below 3 K under 
hydrostatic pressure 0], it can be traced down to 0.1 K 
for Celns-zSnz and an additional first-order phase tran- 
sition TV has been found for 0.25 < x < 0.5 inside the AF 
state |14|. These differences are related to the change of 
the electronic structure induced by Sn doping. Beyond a 
possible tetracritical point at x 0.4 14j, an almost lin- 
ear dependence of Tn{x) is observed. This is in contrast 
to T N (x (x c - a;) 2 / 3 predicted by the 3D-SDW theory 
Q. Thus, the origin of the NFL behavior in this system 
remains an open question, and further thermodynamic 
studies are needed to shed light on the nature of the 
QCP. 

The Celna-^Sn^ single crystals investigated here (0.55 
< x < 0.80) were grown by a Bridgman-type technique. 
Large single crystals with a mass of 15 g were produced, 
analyzed by X-ray powder diffraction and found to be 
of single phase with the proper cubic structure. Within 
the ±2% accuracy of the X-ray diffraction, no impurity 
phases were resolvable. Thin bars with a length 2 < I < 6 
mm, suitable for the dilatometric investigations, were cut 
out. The thermal expansion has been measured in a dilu- 
tion refrigerator using an ultrahigh resolution capacitive 
dilatometer with a maximum sensitivity corresponding 
to Al/l = 10" 11 . 

Fig. 2a shows the volume thermal expansion (3 of 
single-crystalline Celna-^Sn-r with x = 0.55, 0.65, 0.7 and 
0.8 plotted as (3{T)/T vs logT. The volume-expansion 
coefficient (3 is given by (3 = 3 x a, with a being the linear 
thermal-expansion coefficient. For x = 0.55 the broad- 
ened step-like decrease in (3/T at Tn ~ 0.6 K marks the 
AF phase transition, in perfect agreement with specific- 




Figure 2: Volume thermal expansion coefficient f3 of 
Celns-zSn^ single crystals as j3/T vs logT (a). Gray solid 
and black dashed lines indicate T -0 ' 5 and J 1-0 - 4 dependencies, 
respectively. Arrow indicates AF phase transition, (b): De- 
viation of /3/T data for x — 0.65 sample from best power-law 
fits for T < 1 K (squares), T < 2 K (circles), and T < 6 K 
(triangles), respectively, as (/3/T) — (ao + aiT° 2 ) vs logT. 
For clarity, the three data sets have been shifted by different 
amounts vertically. 



heat measurements on the same single crystal |l£| . Upon 
increasing the concentration we find for x = 0.65 and 0.7 
diverging behavior over nearly two decades in T down 
to 80 mK. These data suggest that Tn is suppressed at 
a critical concentration x c « 0.67 ± 0.03, also consistent 
with specific-heat measurements performed on the same 
samples |18| . Finally, for x = 0.8 we recover Fermi- liquid 
behavior, (3{T)/T sa const, for T -> 0. 

In the following, we will analyze the observed NFL 
behavior and make comparison with the predictions of 
the itinerant SDW scenario 0- A best-fit description 
of the x = 0.65 data in the entire temperature range 
0.08 K< T < 6 K according to (3/T = a a + a 1 T a2 reveals 
<i2 = — 0.4±0.01 (see dashed line in Fig. 2a). However, as 
shown in the upper part of Fig. 2b, the deviation between 
the data and this fit shows several broad bumps indicat- 
ing that the fit does not properly describes the data. We 
therefore tried best power-law fits for 0.08 K< T < T max 
with varying T max . For T max — 1 K, the fit is of excel- 
lent quality (cf. solid line in Fig. 2a and lower part of 
Fig. 2b) and the resulting exponent equals —0.5, i.e. the 
value predicted by the 3D SDW scenario jfj. 
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Figure 3: Temperature dependence of the Griineisen parame- 
ter r = Vm/nT • P/C of several Celns-^Sn^ single crystals as 
r(T) vslogT. Vm = 6.25-10 -5 m 3 mol _1 and k t = 1.49-ltr 11 
Pa -1 Hai are the molar volume and isothermal compress- 
ibility, respectively. Arrow indicates AF phase transition. 
The inset displays data for x = 0.65 in a double-logarithmic 
plot. The dotted line indicates the power-law dependence 
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Figure 4: Critical Griineisen ratio r cr = Vm/itr ■ /3 cr /C cr 
for CeIn2.35Sno.65 as logr cr vs logT with critical components 
f3 cr = p(T) - a T and C cr = C(T) - jT derived after sub- 
traction of background contributions (see text). Solid line 
represents F cr oc 1/T e with e = 1.1 ± 0.1. The inset shows 
the low-temperature electrical resistivity of a single crystal of 
similar composition. 



We now turn to the Griineisen parameter defined as 
r = Vm/ Kx ■ (3/C where the constants V m and kt denote 
the molar volume and isothermal compressibility, respec- 
tively. The specific heat has been studied in the temper- 
ature range 40 mK < T < 4 K on the same CeIn 3 _,2;Sn x 
samples used for thermal expansion fl8t Hif. Here, the 
nuclear quadrupole contribution of indium which be- 
comes important below 150 mK has been subtracted. 
Fig. 3 shows a comparison of T(T) for all samples stud- 
ied in thermal expansion. Since the temperature depen- 
dence of specific heat is much weaker compared to that 
of thermal expansion, its influence to the Griineisen pa- 
rameter is rather small. Therefore the variation of T(T) 
for the different Celn 3 _ x Sn x samples is very similar to 
that found in /3/T (compare Fig. 2). Both single crystals 
closest to x c , x — 0.65 and x — 0.7 show a divergent 
behavior down to the lowest accessible temperature with 
very large F values at 0.1 K which are of similar size as 
found for other quantum critical HF systems 

GUIIl- On 

the other hand, saturation is observed for x — 0.55 and 
x = 0.8 being located in the AF ordered and Fermi liq- 
uid regime, respectively. The fact that the divergence of 
r(T) in the quantum critical regime is stronger than log- 
arithmic (compare the double logarithmic representation 
of the x — 0.65 data presented in the inset of Fig. 3) pro- 
vides clear evidence for a well defined (pressure-sensitive) 
QCP in the system. If the disorder present in the system 
would lead to a "smeared" quantum critical regime, T(T) 



could diverge at most logarithmically [9j. 

Another indication for a QCP is the sign change of the 
Griineisen parameter between the ordered and disordered 
regime. As discussed in 22], it is directly related to the 
entropy accumulation near the QCP. The different signs 
of r in the AF and paramagnetic region reflect the oppo- 
site pressure dependencies of the respective characteristic 
energy scales. Below Tn, the effective AF intersite inter- 
action dominates, whose negative pressure dependence 
gives rise to a negative Griineisen parameter, T < 0. On 
the other hand, the positive Griineisen ratio in the para- 
magnetic state is compatible with the positive pressure 
dependence of the Kondo temperature in Ce-based HF 
Systems. 

In order to compare our results for the x — 0.65 sample 
which is located closest to the QCP with the theoretical 
predictions for an itinerant AF QCP 9], we need to cal- 
culate the critical Griineisen ratio F cr (T) oc (3 cr /C cr of 
critical contributions to thermal expansion and specific 
heat. For thermal expansion, /3 cr (T) = j3(T) — a$T with 
ao = 0.3 x 10~ 6 K -2 as determined from the best fit up 
to 1 K, see above. Within the itinerant theory for the 
3D AF case,the critical contribution to specific heat is 
sub-leading 9]: C cr (T) = C{T) - 70T, with C cr < 0, 
C cr /T for T ^ and 70 = C/T | T=0 . For 7o we 
use the value 0.851 Jmol _1 K -2 obtained in 0] from fit- 
ting the low-temperature electronic specific heat in a re- 
stricted temperature range 0.3 K < T < 1.4 K according 
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to C/T = 70 (1 - a'VT). Figure 3 displays a log-log plot 
of r cr (T) versus temperature. We find T cr cx T~ e with 
an exponent e = 1.1 ± 0.1 which is very close to 1, pre- 
dicted by the itinerant theory. Note, that this exponent is 
rather insensitive of 70 subtracted from the specific heat 
data: using 70 = 0.9 Jmol~ 1 K~ 2 and 0.95 Jmol _1 K~ 2 re- 
sults in e = 1.07 and 1.02, respectively. Interestingly, the 
exponent for the critical Griineisen ratio, which theoret- 
ically equals the dimension of the most relevant operator 
that is coupled to pressure @, holds over a much larger 
temperature range than the respective 3D-SDW depen- 
dencies in specific heat and thermal expansion (cf. 
Fig. 2). A similar observation has also been made for 
CeNi 2 Ge 2 Q. 

For those two systems for which an unconventional 
QCP has been proposed, YbRh 2 Si2 and CeCue-zM^ 
(M=Au, Ag), distinctly different temperature depen- 
dences have been observed: T cr oc T~ - 7 in the former 
[lfl | and r cr oc logT [2H in the latter case. It is pro- 
posed in 0, that for magnetically 3D systems without 
frustration the SDW picture should apply. This is con- 
sistent with our Griineisen ratio analysis. 

For the 3D AF case, the itinerant theory predicts an 
asymptotic T 3 / 2 dependence for the temperature depen- 
dent part to the electrical resistivity 0, ■ As discussed 
in 23], the interplay between strongly anisotropic scat- 
tering due to the critical spinfluctuations and isotropic 
impurity scattering can lead at elevated temperature to 
temperature exponents of the resistivity between 1 and 
1.5, depending on the amount of disorder. Systematic 
p(T) studies down to mK temperatures on polycrystalline 
CeIn3_ x Sn x revealed an almost linear temperature de- 
pendence in the quantum critical regime |15| . Similar 
behavior is observed for single crystalline CeIn2.35Sno.65 
as well, see the inset of Figure 4. However, due to the high 
Sn-doping needed to tune the system towards the QCP, 
the resistivity ratio Pzook/Po is of the order of 1 and 
the temperature variation amounts to a few percent of 
Po only, making the comparison with theoretical predic- 
tions very difficult. This indicates that transport exper- 
iments alone are not sufficient to characterize quantum 
criticality in disordered systems. Possibly, also the slope 
of Tn(x) differs from the 3D-SDW prediction because 
disorder is not constant but increases with increasing x. 
However, the algebraic divergence of T(T) for T — > at 
x w x c proves a pressure-sensitive QCP in the system 
and excludes disorder-driven scenarios for the observed 
NFL behavior 0. 

In conclusion, our study on Celns-^Spj; single crys- 
tals by means of the low-temperature thermal expansion 
and Griineisen parameter has proven the applicability of 
the itinerant theory for 3D critical spinfluctuations in 
this cubic system. Since strong contradictions to this 
theory have been found in systems like CeCus.gAuo.i 
or YbRh 2 (Sii_ a; Ge a ;)2 Q with lower crystallographic 
symmetry and, at least in case of the former system, 



strongly anisotropic quantum critical fluctuations, the 
parameter dimensionality obviously plays an important 
role for the nature of HF QCPs. We tentatively classify 
the different HF systems studied by Griineisen analy- 
sis at their res pect ive QCPs as follows: (i) Celna-^Sna, 
and CeNi 2 Ge2 [lfj for which latter system neutron scat- 
tering measurements revealed 3D low-energy magnetic 
fluctuations |24j . show thermodynamic behavior com- 
patible with the 3D itinerant theory, whereas for (ii) 
YbRh 2 (Sii_ ;c Ge x )2 and CeCu 5 . 8 Ag . 2 [U| strong 
contradictions to this model (for both 2D and 3D critical 
spinfluctuations) are observed. Neutron scattering has 
proven 2D quantum critical fluctuations in CeCus.gAuo.i 
[25| while for YbRh2Si2 a complicated behavior with 
competing AF and ferromagnetic quantum critical fluc- 
tuations has been observed |26j . The comparison with 
our results on CeIn3_ 2; Sn a ; suggests that a destruction of 
Kondo screening causing unconventional quantum criti- 
cality is prevented in magnetically 3D systems. 
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